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Abstract 

The variational ansatz for the ground state wavefunctional of QCD is found to capture the 
anti-screening behaviour that contributes the dominant '—4' to the /^-function and leads 
to asymptotic freedom. By considering an SU(N) purely gauge theory in the Hamiltonian 
formalism and choosing the Coulomb gauge, the origins of all screening and anti-screening 
contributions in gluon processes are found in terms of the physical degrees of freedom. The 
overwhelming ant i- screening contribution of '—4' is seen to originate in the renormalisation 
of a Coulomb interaction by a transverse gluon. The lesser screening contribution of '+|' is 
seen to originate in processes involving transverse gluon interactions. It is thus apparent 
how the variational ansatz must be developed to capture the full running of the QCD 
coupling constant. 
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1. Introduction. 



Asymptotic freedom was discovered in some seminal calculations published in the 
1970's. This unexpected result was proved for all non-abelian gauge theories and Quantum 
Chromodynamics (QCD) was born soon after. Khriplovich, calculated the full Green's 
functions of an SU(2) purely gauge theory in the radiation gauge using the spectral 
representation. Within the full Coulomb Green's function of [|T| was written the now well 
known components of the /3— function of the SU(2) Yang-Mills theory. The calculation 
of what is now the QCD /3— function was carried out some years later, and also by 't 
Hooft (unpublished). 



The /3— function for a purely gauge field Yang-Mills theory (no fermions) is, 
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(47r)2 

The contribution '-4' is an overwhelming ant i- screening effect that gives asymptotic free- 
dom. The 'I' is a lesser screening effect. 6*2(6') is the second Casimir operator of the 
SU(N) gauge group in the adjoint representation. 

After the work of and , some attempts were made to develop an intuitive, physical 
understanding of the mechanisms by which the screening and anti-screening effects are 
manifest. In 0, the screening and anti-screening components of the /3— function were 
related to the spin of the gauge field. The anti-screening effect was related to the influence 
of the background field upon the electric dipole density and the screening effect was related 
to the magnetic dipole density. The /?— function for a field of spin S was written. 



'-ll^S r V 
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where the group factor has been omitted. For a gauge field of spin 1, the anti-screening 
effect is —4, as in (|1.1|). 
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The /?— function was decomposed in 0] into the same screening and anti-screening 
components as for a spin 1 particle in [Q, i.e. but with different origins. This 



was achieved from the calculation of the pre-exponential factor, or the renormalisation of 
the charge, for the BPST instanton within the path integral formalism. Expanding the 
action in terms of a deviation from the instanton field, the zero-frequency modes were 
shown to give the anti-screening contribution of '—4' and the positive frequency modes 
the screening contribution of There is no contradiction in the apparently different 
origins of the effects in these calculations since the function can be calculated from 
any physical process and in |^ is related to instantons. 

For a fuller discussion of the calculations and @], the reader is directed to the 
originals or to |^ for a brief overview and discussion in context of the calculation of the 
/3— function of the coupling constant of the variational approach. It is interesting to note 
here, however, that the two calculations, which involve physically different phenomena, 
produce answers for the /3-function that decompose identically into screening and anti- 
screening contibutions. In this paper, the origins of the same screening and ant i- screening 
components will be found in terms of fundamental gluon processes. These gluon process 
are written in the Hamiltonian formalism in terms of the physical degrees of freedom only; 
transverse gluons and Coulomb interactions. 

The corollary of asymptotic freedom is that at low energy there is a strong coupling 
problem, e.g. phenomena such as confinement and chiral symmetry breaking. After a 
quarter of a century of work there is still no complete theory in answer to these ques- 
tions. Despite the suggestion of many promising ideas, the arsenal of necessarily non- 
perturbative methods with which to attack these problems is limited. To have an analytic 
solution for the ground state wavefunctional of an asymptotically free non-abelian gauge 
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theory, with the associated enhanced understanding of the underlying physics, would be 
invaluable for the understanding of these strong coupling phenomena. This was the mo- 
tivation that drove the development of a variational approach to Quantum Field Theory 
(QFT). Although in quantum mechanical problems the variational approach is often easy 
to use - a few qualitative features is usually enough to write an ansatz that will give 
good results for the ground state expectation values - there is rather more complexity to 
overcome to use the approach within QFT, as discussed earlier by Feynman, 0. 

The variational approach was successfully applied to QFT in the exploratory calcula- 
tions of and where QCD and Quantum Electrodynamics in 2 + 1 dimensions (QED3) 
were considered, respectively. Extensive calculations were performed reproducing many 
old and giving new results. The reader is directed to the original papers for details. With 
the confidence given by the many results obtained from the variational approach to QCD 
and QED3, it was necessary to analyze the ansatz for the ground state wavefunctional 
of QCD. In 1^ it was found that the solution of the minimization equation obtained in 
the variational approach results in the variational parameter being fixed away from the 
perturbative value and that the best variational state is characterized by a dynamically 
generated mass scale. The vacuum condensate was also calculated. Hence, the ground 
state wavefunctional captures some known non-perturbative characteristics of the QCD 
vacuum and we would fully expect the /3-function of the coupling constant written in the 
variational ansatz to be like that of QCD. In fact, it was assumed in that it was indeed 
the /^-function of QCD. This assumption was investigated in ^j. 

The calculation of the /3— function from the ansatz for the ground state wavefunctional 
of the variational approach was performed in ^ with the following result. 
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This result is very close to the /3— function of QCD, ( [L . 1|) . The overwhelming anti- 
screening contribution that leads to asymptotic freedom is captured by the variational 
ansatz. The lesser screening contribution is, however, missing. The effective action, ( p.7|) , 
considered in ^ was written in terms of elements of the SU(N) gauge group. Gauss' law, 
which is the only physical constaint imposed upon the ground state wavefunctional gives 
rise to longitudinal gluons with Coulomb interaction and can be used as the generators of 
the group. So, in |Q it was conjectured that the overwhelming anti-screening contribution 
of '—4' is due to the renormalisation of a Coulomb interaction. Further, since all terms 
greater than quadratic were omitted from the Hamiltonian because they gave only small, 
0{g), contributions to the minimization equation, it was conjectured that the screening 
contribution '|' must originate in the interaction of transverse gluons. The aim of the cal- 
culations presented in this paper is to prove this conjecture to be correct and so to obtain 
an intuitive understanding of the origin of the QCD function in terms of fundamental 
interactions in gluon processes and also to show the equivalence of the /3-function of the 
coupling constant of the variational approach to that of QCD. 

In order to prove this conjecture, it is necessary to write QCD in the radiation (also 
known as Coulomb) gauge. In this paper a purely gauge SU(N) Yang-Mills theory in the 
Hamiltonian formalism is fixed in the radiation gauge, diA'} = 0. The time-like and the 
longitudinal components of the gauge field are eliminated leaving only transverse gluons 
and Coulomb interactions in the action. Only physical degrees of freedom appear in the 
Hamiltonian and there are no ghosts. The choice of the radiation gauge, which gives 
transverse gluons and Coulomb interactions explicitly, allows direct comparison with the 
variational approach. 

In section 3, QCD is written in the radiation gauge and the Feynman rules and dia- 
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grams are derived. In the fourth section, the modification of the bare Coulomb interaction 
between two transverse gluons by the insertion of transverse gluon and transvere gluon - 
Coulomb loops, vertex corrections and box diagrams is considered. This is tantamount to 
renormalising the coupling constant of the theory to first order. The /?— function could be 
deduced from any physical process, but to allow direct comparison with the calculation 
of the function of the charge of the variational ansatz we consider a non-local, four 
transverse gluon vertex with an explicit Coulomb interaction, present in the infinite series 
of the action. In the next section, though, details of the calculation of the f3— function 
from the variational approach, [Q, necessary for comparison with QCD are highlighted. 



2. The function of the Variational Approach. 

From the variational ansatz for the ground state wavefunctional of an SU(N) gauge 
theory the /3— function was deduced in ^j. For details the reader is directed to that paper 
and for more understanding of the ansatz for the variational approach to the ground state 
wavefunctional to 0. Here only the necessary details and results will be highlighted to 
enable comparison with QCD in the radiation gauge. 

An SU(N) Yang-Mills gauge theory is described by the Hamiltonian, 



H = d^x 



\Ef + \Bf 



(2.1) 



In the variational calculation of and the analysis of the ansatz in p] all terms higher 
than quadratic in the Hamiltonian were ignored because they give only small contribu- 
tions, of 0{g), to the minimization equation. These terms contained all the information 
about gluon-gluon interactions. 



The ansatz for the ground state wavefunctional, "^[A], was forced to satisfy the con- 
straint of gauge invariance. The gauge transformation Af^ is generated by Gauss' 
law, G''{x)^\J^ = [diE^{x) - gf''^''A''^{x)Ef{x)]<if[A] = 0. Gauss' law gives rise to longitu- 
dinal gluons with a Coulomb interaction. 

In the variational approach the vacuum expectation value of the Hamiltonian is cal- 
culated and minimized. 

<H>= - f (2.2) 



Z 

Since the Hamiltonian is only considered up to quadratic terms and for other reasons 
(calculability, the dominance of a single condensate), the wavefunctional is written as a 
Gaussian. To ensure its invariance under a gauge transformation it is necessary to sum 
over the space of special unitary matrices with the SU(N) group invariant measure. 

^[A] = I DUexp[-^ I d^xd^yAf{x)G-\x - y)Af{y)] (2.3) 

G{x — y) is like a non-local propagator. G is defined to coincide with perturbation theory 
at high momenta and to have a mass gap. The value of the mass gap is the variational 
parameter. 

, fv^ ifP>M2 
^ ^ > \m if P < M2 ^^-^^ 

( |2.3| ) and ( |2.4|) together form the variational ansatz. 

The calculation of the expectation value of an operator (e.g. (|2.2|)) is tantamount 



to performing a three dimensional path integral in Euclidean space. So the exponent of 
can in some sense be considered as an action. Integration over the field A showed 
this action to be a non-local, non-polynomial, non-linear sigma model in three Euclidean 
dimensions, defined by the action r[t/], where, 

T[U] = ^TrlogA^ + ^A"A""A'= 



= Urir^U^diU) 
A'^^(a:,2/) = [G{x - y)5'''' + S''\x)G{x - y)S^"^{y)\~' 
S'^'ix^y) = Ut{t''U+t'U) (2.5) 

are generators of an SU(N) group. Tr is a trace over all indices, tr is a trace over 
colour indices and all summations over indices and integrations over spaces are implicit. 
There are no derivatives in - it makes no contribution to the following analysis and is 
omitted. 

In [Q, this action was analyzed and the /3— function for the coupling constant of the 
theory calculated. The group elements were decomposed into high and low momentum de- 
pendent modes with the ansatz U{x) = Ul{x)Uh{x). Writing Uh{x) = exp[^0"(x)T"] = 
1 + ^(f)"'T'^ — ^{(p^T'^y + 0{g'^) allowed the effective Lagrangian to be written as, 

^L{x,y) = \xiL{x)G~\x-y)Xl^{y) (2.6) 

+ < d,r{x)<P'{y) > G-\x - y)\ldy)] 
-jrj'''>^lLix)G-\x - y)XUy)[< 0^(x)/(x) > 
+ < <l>^{y)Ay) >-2< cj^^ix^'iy) >] 
+^Xl,{x)G-\x - y)xUy)^^^[< > 
+ < 0^(y)0^(y) > -2 < 0^(1/) >]. 

This is an effective low energy Lagrangian. High momentum modes in the region M' < 
k < A were integrated over, where A is the UV cut-off. The scale M' is arbitrary but in 
p[ and here M' = M is chosen for simplicity of calculation and clarity of presentation. 
Hence, 0° depends only upon high momentum modes and Af^^ depends only upon low 
momentum modes. 
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X y 
Fig. 1 

Now the correlations of S"" need to be calculated. 



< (j)^[x) >= (2.7) 

<n^)^\^)>=:^\og[-] (2.8) 

^^^ ^ \ -^^log[M\x-y\] + K M\x-y\<fi ^ ^ 

C and K are finite contributions which are independent of M and A and so are ignored 
in the following. (|2.8|) occurs in in terms such as, 



^Xl,{x)G-\x - y)Xl,{y)^^^ < > . (2.10) 

This can be represented by a Feynman diagram as in Fig. 1, which shows a tadpole 
diagram with a non-local connection. The external lines correspond to low momentum 
fields and the internal solid loop represents the integration over high momentum fields. 
The dotted line corresponds to G~^{x — y) which contains information about transverse 
gluons. In the region k < M, G^^{x — y) = M6{x — y) and the propagator becomes local, 
associating the ends of the dotted line, and the standard tadpole diagram is recovered. 

appears in F^ in terms such as, 

- ^^\i,{x)G-\x - y)\Uy)^^^ < <l>W{y) >, (2.11) 

which can be interpreted in terms of 'horse-shoe' diagram. Fig. 2. The external, internal 
and dotted lines correspond to the same quantities as in Fig. 1, described above. This 
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diagram is unimportant in the region < |x — y| as then it is finite. The interpretation 
of the appearance of |a; — y| must come from the Lagrangian. 

The effective Lagrangian can now be written as, 

T^ix^y) = ^-\l^{x)G-\x-y)\l^{y), (2.12) 



where. 



9 



+ 0(/) (2.13) 



and, 



, M \x - y\ > , , 

I \x-y\ 1-^ y\ ^ M- 



From this the /^-function is easily calculated. 

^) 3 

P{g) = M'^^|,,A = -^AC,{G) + 0{g'). (2.15) 

This result is close to the QCD /5-function calculated in [|l|, 0, and 0] as discussed 
in IP). This prompted the conjecture that the ant i- screening contribution of '—4' to the 
QCD /3-function is due to the renormalisation of a Coulomb interaction. Such Coulomb 
interactions originate with Gauss' law which can be used as the generators of SU(N) gauge 
transformation, i.e. as the generators of U. Further, the conjecture was made that the 
'I' screening contribution must be associated with transverse gluon processes - exactly 
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the interaction terms omitted from the Hamiltonian in the the exploratory variational 
calculations of as they only give small, 0{g), contributions to the minimization equa- 
tion. This conjecture is reliant upon the validity of the original ansatz for the ground 
state wavefunctional. In the next two sections QCD is written in the radiation gauge and 
this conjecture is proved. Thus the equivalence of the /3-function of the coupling constant 
written in the variational ansatz to that of QCD is shown and a physical understanding 
of the origins of screening and anti-screening effects in gluodynamics is obtained. 



3. QCD in the radiation gauge. 

SU(N) Yang-Mills gauge field theory in the Hamiltonian formalism is defined in 
with the constraint of Gauss' law, Diir"- = 0. The covariant derivative is defined as, 

D^w'' = d^w''-gr"'A''^w''. (3.1) 

The generators of the SU(N) internal symmetry, T", and the group structure constants, 
jabc^ have the following properties: 

TrlT^T^] = 

^abcj^abd ^ C'2(G)5'='^. (3.2) 

C2{G) is the second Casimir operator and G denotes the adjoint representation in this 
case. For SU(N), C2{G) = N. Calculations will be performed in Euclidean space. 

Imposing the radiation gauge, diA"; = 0, allows us to write the Hamiltonian density, 
H = + ld,A';d,A^ - gr'^A^A'^A'^^ 
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(3.3) 



From the equivalent to this Hamihonian for the SU(2) Yang-Mills theory in Minkowski 
space Khriplovich ^ deduced the Feynman diagrams and wrote the full Green's functions 
of the theory using the spectral representation in the radiation gauge. This process 
involves decomposing the canonical momenta into transverse and longitudinal components 
and solving to eliminate the longitudinal component. The time-like component of the 
vector field is also eliminated and the Coulomb interaction becomes explicit. The final 
Hamiltonian obtained, with the canonical momenta solved in terms of the components 
of the vector field, is an infinite series in g which contains only transverse gluons and a 
Coulomb interaction. Khriplovich's interpretation was that such an infinite series in the 
coupling constant of the theory was the price for having no unphysical ghost fields. A 
similar procedure will be used for the more complex case of the SU(N) Yang-Mills theory 
in Euclidean space. A diagrammatic treatment of the theory will be used which will 
provide physical insight into the processes from which the screening and ant i- screening 
contributions to the /^-function originate, in contrast to the spectral representation used 
in The Hamiltonian for a purely SU(N) gauge theory shall be written containing only 
transverse gluons and Coulomb interactions. This allows a direct comparison to be made 
between the variational ansatz for the groundstate wavefunctional and QCD. 

To start, the canonical momentum is decomposed into longitudinal and transverse 
parts. Pi and dicj)"' respectively. 



.a 



a 




0. 



(3.4) 
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Thus, we can write the Hamiltonian density, 

1 

4 ■' ■' 1 

Substituting this decomposition of the canonical momenta into Gauss' law gives the 
relation, 

^<^^ = gr'^A%v\ + d,<^^). (3.6) 

This expression can be solved for 0" by a process of reiteration with the definition of the 
inverse Laplacian, 

A-i(a;)/(x) = J d^xA-\x-x')f{x') 
A{x)A'\x-x') = S\x-x'). (3.7) 

0" therefore becomes an infinite series in g, 

r{x) = riix) + r2{x) + rsix) + o{g') 

riix) = gr'^ j d'x'A'\x-x')[A\{x')p'i{x')] 

^\[x) = g^r'T''' J J d^x"A-\x - x')[A\{x')&'^ A-\x' - x")[A''.{x")p]{x")\] 

(j)l[x) = g^f^bcjcdej^efg J ^i^, j ^i^n j (^^^-^(x - x')[A\{x')df ^-'^{x' - x") 

[A%x")df A~\x" - x'")[4(^"'K(^'")]]] (3.8) 
where df = -St- 

By substituting ( |3.8|) into (|3.5|) and identifying p° with —doA"", we obtain the following 
Hamiltonian, 

H = [ - gr'"'diA''A'lA' + — /"''V^'^M^^AMfA"] 
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2 

_9_jabcjade j xd^ x'[A\{x)doA1{x)]A-\x - x')[A%x')doA''^{x')] 

_g3fabcfcdefafg j d'^xd^x' d'^x"[Al{x)d^ A-\x - x')[A'^j{x')doA''.{x')]]A-\x - x") 

[Ai{x")doAi{x")] 

-^g^r^T'^^r^'P'"' j d''xd^x'd''x"d^x"'[A\{x)d'^A-\x - a;')[Aj(x')aoA^"(x')]] 

A-\x - x")[Ai{x")df A-\x" - x"')[A'l{x"')doA'i{x"')]] 
+0(/) (3.9) 

This Hamiltonian contains three terms famihar from the usual Lorentz covariant treat- 
ment of QCD and an infinite series of non-local terms, the non-locality provided by the 
definition of A~^. The transverse gluon propagator, denoted by the conventional 'pig- 
tails', is given by, 

dM = ^^3-^-f- (3-10) 

where — Pi- The superscript of zero denotes that this is the bare propagator. The 
second and third terms give rise to the usual three- and four- gluon interaction vertices but 
in this Hamiltonian these are interactions of transverse gluons only. The three transverse 
gluon interaction with its Feynman rule is shown in Fig. 3. The convention that positive 
momentum flows into a vertex is used and the sum of momenta into a vertex is zero, (i.e. 
in Fig. 3, p -|- g -|- r = 0). 

The fourth, fifth, sixth and all subsequent terms in the action form an infinite series 
in the coupling constant, g. Each term involves a Coulomb interaction between two 
gluons and each successive term differs from the previous by the emission of a transverse 
gluon from the Coulomb interaction. The Feynman diagrams and fully symmetrized rules 
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+igf [(r-q)i^jk+ (q-p)A + (p-r)^ki 



a 1 



Fig. 3 Three transverse gluon interaction 
c^j k 



p+q 



q 



P 



b i 




^1= 



= I f "^f (2p„+q,) (2r,-qo)^,^ ^.A" \q) 



Fig. 4 Four transverse gluon non-local interaction 

corresponding to the fourth, fifth and sixth terms in the action ( p.9|) are given in Figs 4, 
5 and 6, respectively. 



The Coulomb interaction is denoted by a dashed line and is written as, 



ab 



(3.11) 



where, again, = pf and the zero superfix denotes a bare interaction. It is first order 
modifications to the bare Coulomb interaction between two transverse gluons. Fig. 4, 
that are considered in the next section and it is from this that the origins of screening 
and anti-screening effects in gluon processes are deduced. 
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Fig. 5 
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Fig. 6 
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4. Renormalisation of the Coulomb Interaction to First Order 



It is possible to deduce the /3-function from any physical process. The case that shall 
be studied in the following can be regarded as the renormalisation of a four transverse 
gluon non-local interaction. Two transverse gluons meet at a point, x say, and two more 
at another point, x\ with an interaction between the vertices at x and x' . To allow a 
fruitful comparison with the calculation of the /3-function from the variational approach 
a Coulomb interaction between points x and x' is considered. This is distinguishable 
from the case with a transverse gluon exchanged between x and x' by the poles of the 
propagator. Such a four transverse gluon non-local vertex can also be considered as a 
Coulomb interaction between two transverse gluons taking the role of external charges. 
Of course, more complicated interactions involving more Coulomb and transverse gluon 
propagators and vertices are possible. These modifications of the bare propagators and 
vertices are constructed by combining terms of lower order in g (e.g. combining three 
and four transverse gluon vertices with Fig.'s 4 and 5) or from closing loops in terms of 
higher order in g (e.g. Fig. 6). The coupling constant can then be redefined to absorb 
the infinities associated with these modifications and from this the /^-function deduced. 
In this section, all diagrams that contribute to the /3-function up to first order will be 
calculated. 

From previous calculations in QCD, it is known that box diagrams do not contribute 
to the /3-function. Choosing the radiation gauge maps these gluon box diagrams into a 
number of diagrams formed of combinations of transverse gluon and Coulomb interactions. 
Each of these new box diagrams may be individually divergent but summed together they 
cannot contribute to the QCD /3-function. 

The other diagrams that need be considered for the renormalisation of a Coulomb 
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Fig. 7 Tadpole diagram 

interaction to first order are loops in the Coulomb interaction, vertex corrections and 
tadpole contributions. Transverse gluon and transverse gluon - Coulomb loops are both 
divergent. A purely Coulomb loop is not possible by the Feynman rules of the previous 
section. One of the two possible vertex corrections is zero and the other contributes an 
infinity. These cases will be considered in the following subsections. 

First, it is necessary to see that there is no contribution from the tadpole correction 
to the Coulomb interaction. Fig. 7, which is formed from the combination of Fig.'s 3 and 
5. With two ends of Fig. 3 closed to form a loop, the product of a delta function in the 
colour indices and the totally anti-symmetric group structure constant gives identically 
zero. 

4.1. Renormalisation of a Coulomb Interaction by a Transverse Gluon - 
Coulomb Loop. 

The Coulomb interaction can be modified by considering the insertion of a transverse 
gluon - Coulomb loop. This is equivalent to the emission and reabsorbance of a transverse 
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4gQ(G) r(e)i 



Fig. 8 



gluon by the Coulomb interaction. This is exactly the diagram formed by joining the two 
'internal' transverse gluon lines in Fig. 6. 

Fig. 4 shows the four transverse gluon non-local vertex with a Coulomb interaction. 
This diagram is denoted by I. Fig. 8 shows a Coulomb interaction modified by a transverse 
gluon - Coulomb loop between two external transverse gluons. This interaction is written 

as, 

lg'r"'f'"T'''f"^ J ^,'Mn{2p„ + *)(2r„ - ,„)5,j4,A-'(,)A-'(, + s) 
^-'Wf^-is. (4.1) 

where, 



-'ig'C2{G)q^qn^-\q) 



(4.2) 



(2^)4s2(s2 + sg)(s2 + 2sg + g2)- 
It is important to note that since Lorentz invariance has been broken by the choice of gauge 
the time-like component is not treated the same as the spatial components of momentum 
in the loop integration. S is evaluated in Appendix B with the result, 



{Any 



-r(6). 



(4.3) 



where the constants are the same as in ( |2.15| ), r(e) is defined in the standard manner 
and given in Appendix A and t = w — A from the integration with the dimensional 
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5 g^C,(G) r(e)I 
3 (47r)' 



Fig. 9 

regularisation formulae, also given in Appendix A. For QCD in 3 + 1 dimensions e — >^ 0. 
All finite terms are ignored in the calculation of the /9-function. 



4.2. Renormalisation of a Coulomb Interaction by a Transverse Gluon Loop. 



From the combination of two four transverse gluon, non-local vertices. Fig. 4, a box 
diagram and a Coulomb interaction with a transverse gluon loop can be made. The 
fact that box diagrams do not contribute to the /3-function has already been discussed 
above. Fig. 9 shows a Coulomb interaction with a transverse gluon loop inserted. This 
interaction is written as, 

6^''^'d,,y{u) S^'^'dMr + u) _ 
where, 

n = -^C^iG) J 7^(2^0 + r,rA-\r)ND. (4.5) 
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The numerator and denominator, N and D, are, 



N = Si'fSk'i'iu^Sk'j' - Uk'Uf){{r + ufSiqi - (r + w)i'(r + u)i') 



D 



u\u^ + Mo)(r + 2rM + u^){r^ + 2™ + + Tq + 2roMo + Wq) 



(4.6) 



n is evaluated by the lengthy calculations of Appendix C. Throughout the calculation it 
was unobvious that the time-like components of momentum would cancel but there it is 
seen that they do. We finally write the modification to the Coulomb interaction by the 
insertion of a transverse gluon loop as. 



n = — 



3 (47r) 



-r(6) 



(4.7) 



4.3. Calculation of Vertex Corrections 

As well as modifications to the Coulomb interaction by the insertion of loops, cor- 
rections to the gluon - gluon - Coulomb vertex are possible. To first order, two such 
corrections are possible. 

The first correction is obtained by combining the interactions of Fig. 3 and 5 and is 
shown in Fig. 10. This is written as, 

- yr^^/^'V"'7'^^4K2ro - qo)^-\q) J (^(2Po + So)ND, (4.8) 

where. 



N = {-s + q)m5ij'[{s-p)i>Sk'j + {p-q)j5i>k' + {'2p + s-q)k'5i>j] 

{{p + sfdfi^ - {p + s)f{p + s)i/)((s - qfSk'm -{s- q)k'{s - g)^), 

D = [s^(p2 + 2ps + s^)(s^ + 2ps + p'^ + sl + 2poSo + So)(s^ - 2sg + g^) 

{s^ - 2sq + q^ + sl- 2soqo + (4-9) 
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Fig. 11 
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This expression is evaluated in Appendix D. A numerator of 0{s^), where s is the mo- 
mentum variable around the loop, is required for this diagram to be divergent. Although 
this criterion initially seems to be satisfied, all terms of such order cancel exactly. Hence, 
Fig. 10 is convergent and makes no contribution to the QCD /3-function. 

The second possible vertex correction is made by the combination of two triple trans- 
verse gluon vertices. Fig. 3, and a four transverse gluon, non-local vertex. Fig. 4. From 
these components a box diagram is also possible but such diagrams have already been 
discussed above. The second vertex correction is shown in Fig. 11 and written as, 

_ tjaVc'^adejc'cjWgj ^(2.o + go)(2ro - qo)5u^-\q)N D , (4.10) 



where. 



N = Jj/zj// [(2p - s + q)ii5jk" + (-P - 2g - s)fc"5i'j + (2s + g - p)j6ifk"] 
[{2s - p)i5fk' + {-s - p)k'Sij' + {2p - s)j>Sik'] 



[{-p + sf6k"k' - {-p + s)k"{-p + s)k'][s^5jfjii 



[{q + sf5i''i> -{q + s)i"{q + s)i>], 
D = [s^{s'^ + s^)(/ - 2ps + s2)(/ - 2ps + + Po - 2poSo + Sq) 

(g^ + 2qs + s^){q^ + 2qs + s^ + ql + 2goSo + sl)]-\ (4.11) 

This expression is evaluated in Appendix E, where the Coulomb interaction between two 
transverse gluons with a transverse gluon correction about one of the gluon - gluon - 
Coulomb vertices is finally written as IF, with. 

Two such vertex corrections are possible to first order, one at either end of the Coulomb 
interaction. So, the modification to the Coulomb interaction between a pair of external 
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transverse gluons to second order in g is, 



/(s + n + 2r) = J 



(4.)^ 3 



(4.13) 



With some normalisation, a redefinition of tlie coupling constant written in I is seen 
to lead exactly to the QCD /5-function. The anti-screening contribution of '—4' is seen 
to come from the renormalisation of a Coulomb interaction by a transverse gluon. This 
is exactly the sum of Fig.'s 1 and 2 from the variational approach, where contains 
information about transverse gluons. The contribution is seen to arise from a screening 
contribution of '|' from a transverse gluon loop and an anti-screening contribution of ' — |' 
from two vertex corrections. These contributions come from the interaction of gluons, 
exactly the terms omitted from the exploratory variational calculation of [0. 



5. Conclusion 

In the exploratory variational calculation of [0, the assumption was made that the 
coupling constant written in the variational ansatz for the ground state wavefunctional 
of QCD ran as the coupling constant of QCD. This assumption was found to be nearly 
correct in ||^ where the /3-function of coupling constant of the variational ansatz was cal- 
culated. It was found to have only a '—4' ant i- screening contribution and not both the 
known ' — (4 — |)' screening and anti-screening contributions of QCD. The only physical 
constraint imposed upon the variational ansatz was Gauss' law. The /3-function from the 
variational approach was calculated from a non-linear, non-local sigma model in three Eu- 
clidean dimensions written in terms of elements of the gauge group for which Gauss' law 
can be used as the generators. Gauss' law gives rise to longitudinal gluons with Coulomb 
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interaction. Hence, the conjecture was made in ||^ that the anti-screening '—4' contribu- 
tion must come from the renormahsation of a Coulomb interaction and the screening 'i' 
contribution must arise from gluon interactions - exactly the non-linear terms missed out 
of the original calculations, 0. 

This conjecture, however, rested heavily upon the validity of the original ansatz. In 
order to prove this conjecture to be correct, it was explicitly checked in this paper by 
performing calculations for QCD in the Hamiltomian formalism and choosing the radiation 
(or Coulomb) gauge. The verification of the conjecture suggests that the original ansatz 
captures the anti-screening behaviour of gluons in QCD that leads to asymptotic freedom 
and is essential to address strong coupling problems. Further, a full understanding of 
the origins of screening and ant i- screening effects in terms of gluon interactions has been 
obtained. 

The anti-screening contribution of '—4' is seen to come from the renormahsation of a 
Coulomb interaction by a transverse gluon. This is shown to be exactly equivalent to the 
anti-screening contribution obtained from the variational approach. The 'I' contribution 
is seen to arise from a screening contribution of '|' from a transverse gluon loop and 
an anti-screening contribution of ' — |' from two vertex corrections. These contributions 
come from the interaction of gluons, exactly the terms omitted from the exploratory 
variational calculation of 0. It is hoped, that by a more refined consideration of the 
initial ansatz to include the non-linear terms that correspond to gluon interactions, it can 
be written to capture the exact running of the QCD coupling constant and can be used 
as a non-perturbative tool to investigate the strong coupling phenomena. 
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Appendix 



A. 



The Gamma function and dimensional regularisation formulae and their properties are 
widely published, e.g. [l^, and so shall not be given here in detail. The dimensional 
regularisation formulae are generally only quoted up to a numerator of where / is 

the momentum integration variable. In Appendices D and E the formulae for numerators 
of 0{l^) and 0{l^) are required. The necessary formulae are given below. 



J2wl 111111 

+Sup{S^aSal3 + ^aa^pp) + ^pai^pa^u^ + ^ua^pji) 

V[A — ?) — w) 

d^^l lJpL,kfl^ 1 63 ^ , ^ ^ ^ ^ ^ ^T{A-A-w) 



(^""1 IJpPPP _ 1 315 r{A-A-w) 

(2^)2- (/2 + M2)^ ~ (47r)"'r(A)T" "'^ (M2)^-4— ^ ■ ^ 
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B. 



In this Appendix, E of Section 4.1 is calculated. 

The loop integration over momentum requires the introduction of three Feynman param- 
eters, one of which is eliminated to give, 

J (27r)4 s^{s^ + sl){s^ + 2sq + q^) ~ Jo Jo ^ ' J (27r)4 + sgxs + M^f ' 

(B.2) 

So = So 

s = s + g(l - xi - a;2) 

= g2^i_a;^_a;2)(a:i + X2) (B.3) 

Using the standard dimensional regularisation formulae of Appendix A, integration over 
the spatial and time-like components separtely gives, 

-r(6), (B.4) 



VC2(G'). 



(47r)^ 



c. 



In this appendix, 11 of Section 4.2 is evaluated, where, 

n = -^^C,{G) J ^(2^0 + rorA-\r)ND. (C.l) 
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The numerator and denominator, N and D, are, 

N = Si'fSk'i'iu^Sk'j' - Uk'Uf){{r + ufSiqi - (r + u)ii{r + u)ii) 



D 



u\u^ + Mo)(r + 2ru + u^){r + 2™ + m'' + Tq + 2roMo + Wq) (C-2) 



-1 



The method employed in the evaluation of H is similar to that used for the evaluation 
of S in section 3.1 although more complex. It is necessary to introduce four Feynman 
parameters to use the methods of dimensional regularisation. Doing this and performing 
a shift in momentum, the denominator can be written, 

_ fl fl—xi rl—xi—X2 

D ^ dxi dx2 / dx3rU)\u'^ + ul(l - xi - X3) + M^]"^ 

Jo Jo Jo 

u — u + r{l — xi — X2) 

ro(l - xi - X2 - X3) 
Uq = uo-{ 

(1 - Xi - X3) 

= r^(l — xi — X2) + ^0(1 — xi — X2 — X3) — r^(l — Xi — ^2)^ 
rg(l -X1-X2- XaY 

(1 - xi - X3) ■ ^ ■ ^ 

The momentum in the numerator must also be shifted but before this is written it is useful 

to consider what powers of uq and Ui in the numerator will contribute to the /3-function. 

Prom the formulae of the Appendix we note that if the numerator is quadratic in uo then 

the denominator must be at most of power |. i.e. 



3 -r(e) (C.4) 



In order to have this second integration, because the initial integrand has a denominator 
of power 4 then the numerator must have a component u^. i.e. 



/ 



(C.5) 



(27r)3 [u2 + ^2 ^ ^2] f [^2 ^ ^2] f 

Similarly, if there are no powers of uq in the numerator, is required. All odd powers 
of Uq and Ui integrate to zero. Therefore, only terms quadratic and quartic in Ui need be 
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considered in the shifted numerator, N. 

N = 2u^u^ + u^{r — 2fY + UjUi{—Qfj{r — f)i + (r — f)j{r — f)i + fjfi) + ... 
The modification to the Coulomb is now written as, 

/■I /■! — xi r-l—x\—X2 

/ dxi / dx2 / dx^r{A)N[u^ + ul{l - Xi - Xg) + M^]" 

= n^i + + (r - 2f)'(nBi + n^s) 

+ (-6fj(r - f)i + (r - f)j(r - f)/ + fjfi){Il.ci + nc2), 
where II^i to Ilc2 will be evaluated separtely. 

^ r(-i) 

(4^)2r(4)(l-xi -a;3)t 

^A2 = 2{ro-—^ f f^u^u^[u^ + ul(l-:c^-xs) + M']- 
{l-Xl- X3) J {27rp 

15(ro - TT^^y 

_ V " (1-3:1-0:3)^ p/ N 

2(47r)2r(4)(l-xi -X3)5 ^' 
/ + ^^^^ " " + ^'1"' 

^ r(e), 



n 



Bl 



{Anyr{A){l - X, - xsf^ 

(l-xi-xa)^ i (27r)4 



finite 

d'^u 
(2^ 



/d 11 - 
-7—riAulu,ui[u^ + ul{l - xi - X3) + M^]-^ 



(47r)2r(4)(l-Xi-X3)^ 
= finite. 
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Using the definition of tlie Gamma function in the Appendix, r(e) = — r(— 1 + e). All 
finite terms are ignored. ( |C.7| ) is now written as, 

15M2 



n 



C2(G) _i ri-xi i'i-xi-x2 

■ — j^;^T{e)A (r) j dxi / dx2 / rfxg 



15/ 2ro 

2 V° (i-xi-xa) 



JO 
2 

(1 - Xi - Xs) 







+ (1 — Xi — X3) 2 (^4(r — fY — 12f(r — f) + 4rr 



(C.9) 

Where spatial indices are omitted, summation is implied. Considering terms separately, 
again, execution of the Feynman parameter integrals gives the results below. 



1 rl-xi j-l-xi-X2 — 15Af^ 

dxi j dx2 / dx'i 3- 

"'0 JO (1 - Xi - X3)"2 



1 rl~xi i'l-xi-X2 X5 

dxi / dx2 / dx3 — 

Jo Jo 2 



2fo 



26 2 2 2 
1 2 



;i - xi - X3) 



'1 - xi - X3) 2 = -r" 



l—x-i 



1 /•1-XI-X2 3 

(iXi / (iX2 / (ix3(l — Xi — X3)~2 

Jo Jo 



4(r — r) — 12r(r — r) + 4rr 
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-r 



(C.IO) 



We finally write the modification to the Coulomb interaction by the insertion of a trans- 
verse gluon loop as, 

(C.ll) 



3 (47r)2 



D. 

In this section one of the vertex corrections of Section 4.3 is evaluated and found to be 
convergent. Fig. 10 is represented by the expression, 

- yr^^/^''r'7'^%K2ro - qo)A-\q) J ^(2po + So)ND, (D.l) 
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where, 

N = {-s + q)mSij' [{s - p)i>5k'j + {p- q)jSi>k' + {2p+ s- q)k'Si'j] 

{{p + sf5j,i> -{p + s)ji{p + s)i>){{s - qfdk'm -{s- q)k'{s - g)m), 

D = [s^{p'^ + 2ps + s'^){s'^ + 2ps + p'^ + sl + 2poSo + sl){s'^-2sq + q'^) 

(^2 _ 2sq + q^ + sl- 2soqo + qoT' ■ (D-2) 

Following the method of the previous sections, this denominator requires five Feynman 
parameters to allow the use of the dimensional regularisation formulae. It can be written, 
after a shift in momentum, as, 

fl rl—xi — X2—X3 



D = / dxi... dx4T{5)[s'^ + sl{l-xi~X2-X4) + M'^]-\ 

Jo Jo 

s = s+p, 

, Po 

So = So + -, 

(1 - Xi - X2 - X4) 

M2 = p2(^a,2 + X3) + q'^{l - xi - X2 - X3) + pIxs + 5o(l - Xi - X2 - X3 - X4) 

[ \ — Xi — X2 — X4) 



Inserting ( |D.3| ) into (p.2| ) and performing the same shift in momentum throughout the 
integrand, in a similar manner to that used in section 3.2 around ( C.4 ) and ( |C.5|) it is 



possible to evaluate what powers of u and uq in the numerator contribute to the /^-function. 
First, we see that the term linear in sq integrates to zero by symmetry. Secondly, for the 
final integration over time-like momentum to be divergent with a constant numerator, a 
denominator with maximum power | is required. For an initial denominator of power 5 
this implies a numerator with terms of O(s^). Writing out the terms in N, 

N = -Sm[Si'Sk'j + Sk'Si'j][s'^Sii' - SiSi'][s'^Sk'm - Sk'Sm] + 0{s^) 

= 0{s^) (D.4) 
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Therefore, Fig. 11 makes no contribution to the QCD /3-function. 



E. 



In this section the second vertex correction of Section 4.3 is calculated. This corresponds 
to the evaluation of Fig. 11, which is represented by the expression, 

- ^lr^''^'r'^r''^9fW, J ^(2^0 + go)(2ro - qo)5,iA-\q)ND, (E.l) 



where. 



N = Sf'i"[{2p - s + q)i'Sjk" + {-p -2q- s)k"Si>j + {2s + q - p)j5i'k"] 

[{2s - p)iSj>k' + {-s - p)k'Sij' + {2p - s)f6ik'] 

[{-p + sf5k"k' - {-P + s)k"{-p + s)k'][s'^Sfj" - SfSj"] 

[{q + s)^5i/v - {q + s)in{q + s)i>], 
D = + si) {/ - 2ps + s^) {p^ - 2ps + + Po - 2poSo + Sq) 

{q^ + 2qs + s^){q^ + 2qs + s^ + q^ + 2qoSo + sl)]-\ (E.2) 

This denominator requires six Feynman parameters to allow the use of the dimensional 
regularisation formulae. It may be written as, 

D = / dxi... dxr,r{6)[s'^ + §1{1 - xi - X3 - X5) + M^]-^ 

Jo Jo 

s — s — p 

Po 



So = So 



{1 - Xi - X3 - X5) 

p = p{xi + X2) - q{l - xi - X2 - X3 - X4) 
Po = P0X2 - qo{'^ - xi - X2 - X3 - X4 - X5) 
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= + X2) + pIx2 + q^{l - xi - X2 - - X4) 

+go(l - xi- X2- X'i- xa - xs) -p^ - 



(E.3) 



Xi - 2:3 - Ssj 

Performing the same shift of momentum in the rest of the integrand, (Pi).l|) can be written 
as, 



y jab' c' jade jc'cg jbb'g 



JO 



1 — Xi — X2 — X3 — 2:4 



dx^ 



dH 



■(2ro - go)4« 



2poa:2 - go(l - a^i - 2x2 - 0:3 - 2x4 - xs) 

^■SO H r 

(1 - Si - 2:3 - Xs) 

[g2 + 5^(1- xi-xa- 0:5) + MT'- 



(E.4) 



The term hnear in sq integrates to zero by symmetry. Using methods described above, 
all the remaining terms in the numerator are constant in sq and so only terms of 0{s^) 
will be infinite. The momentum shifted numerator can be written as, 



N = %s,~sfs° + 0{s'). 



(E.5) 



Within 

A) is, 

8 



the integration over the loop momentum (using the formulae of Appendix 



d^S 



S i S j S S S 



315- 



(47r)2(l - xi - X3 - X5) 



(27r)4 [§2 + 52^1 _ _ 2:3 - X5) + M2]6 
The Feynman parameter integrals give, 

2poa^2 - go(l - xi- 2x2 - 3^3 - 2x4 - 3:5) 



rr(e). (E.6) 



1 r\~x\—X2 — xz—X4, 

dxi / dx^ 
Jo 



;i - xi - X3 - X5) 



^(2po + go) 

(E.7) 



Therefore, substituting ( Pi).6|) and ( [E.7| ) in (|E.4| ), the vertex correction can be written as, 
2 9' 



jabcjadejccgjbbg^^p^ + go)(2ro - go)%4^ (g)! (e) 



(E.8) 



3(47r)2- 

Now, the four structure constants need to be shown to be equal to a product of two. As 
defined in (p.2|), and /"'"^ are the generators and structure constants of an SU(N) Lie 
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group, respectively. In the adjoint representation, 

This allows relults known, e.g. 0, for the generators of an SU(N) Lie group to be 
exploited. For instance, the result Tr(rT''T^) = yz/"'"^ allows the evaluation, 

jayc'jadejc'cgjWg ^ _,(T«),,,, (r'=),,,(T^),,,/'^'^^ = f-'^ f^'^^ (E.IO) 

Therefore, the Coulomb interaction between two transverse gluons with a transverse gluon 
correction about one of the gluon - gluon - Coulomb vertices is finally written as IF where, 

2 g^C,{G) 
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